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Abstract 


Introduction 


A computational and experimental study was 
conducted to investigate the effects of multiple 
injection ports in a two-dimensional, convergent- 
divergent nozzle, for fluidic thrust vectoring. The 
concept of multiple injection ports was conceived to 
enhance the thrust vectoring capability of a convergent- 
divergent nozzle over that of a single injection port 
without increasing the secondary mass flow rate 
requirements. 

The experimental study was conducted at static 
conditions in the Jet Exit Test Facility of the 16-Foot 
Transonic Tunnel Complex at NASA Langley Research 
Center. Internal nozzle performance was obtained at 
nozzle pressure ratios up to 10 with secondary nozzle 
pressure ratios up to 1 for five configurations. 

The computational study was conducted using the 
Reynolds Averaged Navier-Stokes computational fluid 
dynamics code PAB3D with two-equation turbulence 
closure and linear Reynolds stress modeling. Internal 
nozzle performance was predicted for nozzle pressure 
ratios up to 10 with a secondary nozzle pressure ratio of 
0.7 for two configurations. 

Results from the experimental study indicate a 
benefit to multiple injection ports in a convergent- 
divergent nozzle. In general, increasing the number of 
injection ports from one to two increased the pitch 
thrust vectoring capability without any thrust 
performance penalties at nozzle pressure ratios less than 
4 with high secondary pressure ratios. 

Results from the computational study are in 
excellent agreement with experimental results and 
validates PAB3D as a tool for predicting internal nozzle 
performance of a two dimensional, convergent- 
divergent nozzle with multiple injection ports. 
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A desirable goal of a fighter aircraft designer is to 
increase the agility, maneuverability, and survivability 
of the jet aircraft. Thrust vectoring dramatically 
increases the agility, maneuverability, and survivability 
of jet fighters 1 ' 7 . For example, augmentation of 
conventional controls with pitch and yaw vectoring was 
studied and showed to allow operation in the post- stall 
regime, where conventional controls are ineffective 1 ' 2 . 

Thrust vectoring on fighter aircraft is currently 
achieved mechanically with movable flaps, which 
directs flow exiting the nozzle. The F-15 SMTD, F-18 
HARV, and F-22 are all examples of fighter aircraft 
that employ mechanical thrust vectoring 8 ' 10 . There are 
several problems with movable flaps for thrust 
vectoring. Movable flaps require heavy mechanical 
hardware, which adds weight to the fighter and 
increases nozzle complexity, adding to aircraft 
maintenance requirements. In addition, movable flaps 
can increase the signature of the nozzle, thus reducing 
the survivability of the aircraft. 

An alternative to mechanical thrust vectoring is 
shock vector fluidic thrust vectoring. The shock vector 
fluidic thrust vectoring concept is shown in figure 1 . 



Figure 1. Shock vector fluidic thrust vectoring. 


Shock vector is achieved by fluidic injection into the 
divergent portion of the nozzle, which causes a shock in 
the primary flow, and turns the flow supersonically. 
Since there is no mechanical hardware other than 
control valves associated with fluidic thrust vectoring, 
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the problems associated with movable flaps are 
eliminated. However, fluidic thrust vectoring 
introduces new problems. The main problem is that the 
secondary stream draws air from the primary air supply, 
which reduces the maximum possible thrust that can be 
achieved by the engine. If the amount of secondary air 
drawn from the primary air supply is reduced to 
maximize the thrust of the engine, then the fluidic thrust 
vectoring may be minimized, and the benefits of fluidic 
thrust vectoring may be negligible when compared to 
movable flaps. The concept of utilizing multiple 
streams for fluidic thrust vectoring was conceived to 
maximize the pitch thrust vector angle without 
increasing the secondary flow requirement. 

The goal of this study was to investigate the effects 
of multiple injection streams for fluidic thrust vectoring 
in a two dimensional, convergent-divergent (2DCD) 
nozzle. In addition, a comparison between 

experimental and computational results was made to 
validate PAB3D as a viable tool for predicting nozzle 
flows with multiple injection streams. This study was 
not meant to produce optimal configurations for 
multiple injection streams for fluidic thrust vectoring, 
but is meant to show the viability of multiple streams 
for fluidic thrust vectoring. 

Five experimental configurations with secondary 
flow were tested. Internal nozzle performance was 
obtained at nozzle pressure ratios (NPR) up to 10 with 
secondary pressure ratios (SPR) up to 1 for all five 
configurations, which corresponds to secondary weight 
flow up to 6% of the primary mass flow rate. A 
computational grid of two experimental configurations 
was developed and internal nozzle performance was 
predicted for NPR up to 10 with an SPR of 0.7 (4% of 
primary mass flow rate) for the two configurations. 


Nomenclature 

F a axial force, lbf 

Fj primary flow ideal isentropic thrust, lbf 

F n normal force, lbf 

F r resultant force, lbf 

F R /Fj resultant thrust ratio 

k turbulent kinetic energy, Btu 

p fluid static pressure, psi 

p a ambient pressure, psi 

Pt, P ort total pressure of secondary flow, psi 

p t j total jet pressure, psi 

x axial direction, in 

x t axial location of throat, in 

y normal direction, in 

y + law-of-the wall coordinate 

6 p pitch thrust vector angle, deg 

8 turbulent energy dissipation 



Abbreviations 

2D 

two dimensional 

config 

configuration 

taj 

injection 

CD 

convergent-divergent 

CFD 

computational fluid dynamics 

EXP 

experimental data 

NPR 

nozzle pressure ratio 

npr d 

design nozzle pressure ratio 

SPR 

secondary pressure ratio 
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Experimental Arrangement 
Test Facility 

The experimental study was conducted at static 
conditions in the Jet Exit Test Facility of the 16-Foot 
Transonic Tunnel Complex at NASA Fangley Research 
Center. Models are mounted to the dual-flow 
propulsion simulation system on a test stand located in 
the test bay of the facility. The model is supplied with 
a regulated continuous flow of clean, dry, and stable air. 
The nozzle flow is exhausted to atmospheric pressure 
inside the test bay and then vented through two 
mufflers on the roof. The test is monitored by a closed- 
circuit television in an adjacent data acquisition room. 

Dual-Flow Propulsion Simulation System 

A schematic of the dual-flow propulsion system is 
shown in figure 2. 



The dual-flow propulsion simulation system is 
connected to two separate heated air supply systems 11 
with flow that can be controlled at constant total 
temperatures between 60 °F and 95 °F and weight flow 
rates between 0.1 lb f /sec and 20 lb f /sec. One supply is 
used for primary flow and the other for secondary flow. 
Both of the air systems were utilized for this 
experiment. The propulsion system uses plenums, 
choke points, and turns to reduce the effects of 
momentum transfer and pressurization as air is 
delivered to the nozzle. 

Test Nozzle 

A photo of the test nozzle is shown in figure 3. 
The nozzle used in this study was a non- axi symmetric, 
2DCD nozzle with a design nozzle pressure ratio 
(NPR d ) of 8.78, a nominal throat area of 4.317 in 2 , and 
a width of 3.990 in. The upper divergent flap hardware 


of the nozzle was designed to accommodate fluidic 
injection. In addition, the nozzle was equipped with 
full length, optical quality plexiglass sidewalls for 
shadowgraph flow visualization as shown in figure 3. 
Details of the nozzle flap contour are defined in figure 
4. 



Figure 3. Test nozzle. 
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X 
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0.000 

0.000 

B 

0.000 

■0.614 

C 

0.000 

1.386 

D 

0.917 

1.163 

E 

1.988 

0.611 

F 

2.394 

0.553 

G 

2.430 

0.559 

H 

2.275 

1.166 

1 

4.550 

0.972 


Figure 4. Nozzle flap contour. 

The test nozzle was designed with interchangeable 
divergent flaps in order to investigate the effects of 
varying the distance between two injection ports in the 
upper divergent section. A total of five multi-port 
injection configurations were tested. The divergent flap 
geometry with multi-port injection and a table of 
important port geometry parameters for each 
configuration are shown in figure 5. 

Instrumentation 

Force and moment data on the model were 
measured using a six-component strain-gauge force 
balance located on the dual-flow propulsion simulation 
system (see figure 2). The flow rate of each air system 
was measured with a multiple critical venturi located 
upstream of the dual-flow propulsion simulation 
system. 
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Configurat ion 

Upst ream Slot 

x, in 

lengt h, in 

widt h, in 

area, in z 

1 

4.10 

3.49 

0.08 

0.28 

2 

3.87 

3.49 

0.04 

0.14 

3 

3.64 

3.49 

0.04 

0.14 

4 

3.41 

3.49 

0.04 

0.14 

5 

3.19 

3.49 

0.04 

0.14 


Note: The Downstream Slot geometry is constant for all five configurations; 
x = 4.10 in, length = 3.49 in, width = .04 in, and area = .14 in?. Also, 
the axial position (x) of the slots is measured from the centerline of the 
slot. 


Figure 5. Divergent flap geometry. 


The total or reservoir conditions of the nozzle were 
measured in the nozzle instrumentation section (see 
figure 2). The upper convergent flaps of the nozzle 
were equipped with 7 static pressure taps at the 
centerline. The upper and lower divergent flaps were 
equipped with up to 82 static pressure taps located at 
the centerline and at 0.40 inch from the nozzle sidewall. 
The secondary plenum on the model was equipped with 
1 total pressure probe at the centerline and 4 static 
pressure taps located 0.5 inch and 1.25 inch from either 
side of the centerline. All pressures probes and taps 
were connected to electronically scanning pressure 
modules. 


Some experimental results are shown with error band to 
compare the data with computational fluid dynamics 
(CFD) predicted data. Details of the error band 
calculation for the experiment can be found in reference 
13. 

Flow Visualization 

A shadowgraph system was used to visualize the 
shock patterns of the flows in this study. The system 
utilizes a point light source, which is transmitted 
through the test section (plexiglass sidewalls of nozzle). 
The emerging light is refracted due to the flow field in 
the test section and captured in a vertical plane placed 
at a distance behind the test section, creating a shadow 
effect 14 . 

Paint-oil flow visualization was used to visualize 
flow patterns along the nozzle surfaces. The paint-oil 
flow mixture was made with a combination of oil paint 
and a few drops of linseed oil to modify the stiff oil 
paint into a thick consistency that flowed smoothly. A 
thick paint consistency was required to avoid complete 
removal of the paint under supersonic flow conditions 
in the divergent section of the nozzle. 


Computational Fluid Dynamics Study 


Data Acquisition and Reduction 

Data from the test nozzle, propulsion simulation 
system, and air supply system were recorded 
simultaneously. To achieve steady-state data, 50 
frames of instantaneous data were taken at a sampling 
rate of 10 Hz. The measurements were averaged and 
corrected using the data reduction equations found in 
reference 12. 

The following conventions were used for data 
reduction. Nozzle pressure ratio (NPR) is the ratio of 
jet total pressure p tj to atmospheric pressure p a . 
Secondary pressure ratio (SPR) is the ratio of the total 
pressure of secondary flow p tport to the jet total pressure 
p t j . Resultant thrust ratio (F R /Fj) is determined from the 
ratio of the resultant gross thrust F R to the ideal, 
isentropic thrust F { (calculated from the measured 
weight flow and stagnation conditions using one- 
dimensional theory). The resultant gross thrust F R is 
determined from the square root of the sum of the 
squares of measured axial F A and normal F N forces. 
Lastly, pitch thrust vector angle 6 p is the angle 
calculated in the pitch thrust-vector plane at which the 
resultant gross thrust vector is deflected from the nozzle 
axis and is determined from equation 1 . 

d p = tan” 1 L- (1) 

^ A 


The Reynolds-averaged Navier Stokes CFD code 
PAB3D was used with two-equation k-e turbulence 
closure and linear Reynolds stress models to predict the 
internal performance of configurations 1 and 5 from the 
experimental study (see figure 5). 

Flow Solver 

PAB3D solves the simplified, Reynolds-averaged, 
Navier-Stokes equations in conservative form, obtained 
by neglecting streamwise derivatives of the viscous 
terms. PAB3D has three viscous models: coupled, 

uncoupled, and thin layer Navier Stokes models. For 
this study, the thin layer model was used. The Navier- 
Stokes equations are solved with one or more of the 
following numerical schemes: the flux vector- splitting 
scheme of van Leer 15 , the flux difference-splitting 
scheme of Roe 16 , and a modified Roe scheme. 
Typically, Roe's upwind scheme is used to evaluate the 
explicit parts of the Navier-Stokes equations and van 
Leer's scheme is used for the implicit part. The inviscid 
terms of the Navier-Stokes equations are split and 
upwind differenced, while the diffusion terms are 
centrally differenced. The details and applications of 
these methods are given in references 15 to 17. For the 
present 2D study, it was possible to solve the entire 
problem implicitly, so van Leer's scheme was utilized, 
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which speeds convergence and reduces computational 
time. 

The PAB3D solver is fully implicit in the cross- 
stream direction. Normally, three-dimensional 
solutions are developed with a relaxation procedure that 
sweeps downstream and upstream with an implicit 
procedure to update each cross-stream plane. In a 2D 
domain, there is only one cell in the cross-stream plane. 
Therefore, this relaxation sweep is not required for a 
two-dimensional computational domain because the 
streamwise plane can be swapped with the cross- stream 
plane to obtain a fully implicit domain. The index 
swapping procedure typically increases the rate of 
convergence and decreases the computational space and 
time required. 

Turbulence Model 

Turbulence modeling was required to predict 
solutions for the flow field of this study. PAB3D can 
perform several turbulence simulations by 
implementing either a 2-equation, linear or nonlinear 
model. A 2-equation k-8 turbulence model was used 
for this study based on previous work done on a 2DCD 
nozzle 18 . 

The k-8 model has a singularity at solid surfaces, 
so either a damping function or a wall function must be 
used to capture the turbulent behavior at solid surfaces. 
Since a wall function is not appropriate for separated 
flows, a modified Jones and Launder form 19 of the 
damping function was utilized at the wall because 
separated flow regions along the divergent section were 
expected at over-expanded conditions. A high 
Reynold's number model of Jones and Launder with no 
damping function was implemented in the free stream 
blocks. 

Computational Domain 

Figures 6(a-c) and 7(a-b) show the computational 
domains for configurations 1 and 5, respectively. A 
total of 364812 and 461100 grid cells defined the 
computational domain for configurations 1 and 5, 
respectively. Relative to the nozzle exit, the ambient 
region surrounding the nozzle extended approximately 
19 throat heights downstream, 23 throat heights 
upstream, and 126 throat heights normal to the jet axis. 
The first cell height along the inside of the nozzle was 
approximately y + =0.5 and there were approximately 40 
cells in the boundary layer grid for each configuration. 



(a) Far -field and nozzle. 



(b) Close-up of nozzle. 


(c) Close up of injection port with primary flow block. 
Figure 6. Computational domain for configuration L 
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Initial and Boundary Conditions 


A total of seven boundary and initial conditions 
were used to properly constrain and initialize the flow. 
The initial conditions consisted of a Mach number 
applied to the inflow duct of the nozzle and to the 
secondary flow plenum to initialize the flow. A wall 
"trip" point, located near the beginning of the inflow 
duct, was used to initialize the turbulent boundary layer 
of the flow. The boundary conditions are shown in 
figure 8. A stagnation condition was applied to the 
inflow duct of the nozzle and to the plenum of the 
secondary flow. The stagnation condition was chosen 
to match experimental conditions for total temperature 
and pressure for the inflow duct and secondary plenum. 
The static ambient region surrounding the nozzle was 
defined by various conditions. First, the left face was 
defined by a subsonic inflow condition. On the top and 
bottom faces, a characteristic boundary condition was 
defined. Lastly, a smart boundary condition that 


switched between constant pressure outflow (subsonic) 
and first order extrapolation (supersonic), depending on 
the local Mach number, was applied to the right face. 
All solid walls were treated as no-slip adiabatic 
surfaces. 


Characteristic BC 


4 Subsonic Inflow 

Stagnation 
\ \ 


y_i\ 

Smart BC 

1 L<LU 

'""f T 

No-slip Adiabatic 
Wall 



Characteristic BC 
Figure 8. Boundary conditions. 


Nozzle Performance Calculation 

Two performance codes were utilized, POST 20 and 
Oracle 21 , to calculate aerodynamic quantities for the 
computational study. POST was used to generate the 
raw data for the nozzle such as Mach number and 
pressure. Oracle was used to define the control volume 
used to calculate the aerodynamic forces and propulsion 
coefficients. 

The performance codes were used to generate 
internal surface pressure data, thrust ratio, and pitch 
thrust vector angle data. In addition, the codes were 
used to construct Schlieren flow visualization images of 
the nozzle flow. Schlieren flow visualization is very 
similar to shadowgraph flow visualization and was used 
to compare with the experimental shadowgraph images. 


Experimental Results 

Experimental data was taken on five configurations 
(see figure 5) at NPR ranges from 2 to 10 with SPRs of 
0, 0.4, 0.7 and 1.0. Internal flow features (static 
pressures and flow visualization), pitch thrust vector 
angles, and thrust performance will be presented here. 
The discharge coefficient will not be presented because 
the plexiglass sidewalls of the nozzle flexed slightly 
under pressurization and an accurate measurement was 
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not possible. Reference 13 has a detailed discussion on 
the discharge coefficient measured in this experiment. 

Internal Flow 

Internal flow features are shown for configuration 
1 (baseline injection) in figures 9 to 13. Figure 9 shows 
the centerline, upper surface, static pressures for 
configuration 1, SPR=0.7 (4% of primary mass flow 
rate). 



NPR 

— 2.01 95 
HU- 2.5008 
- 9 - 3.01 05 
3.5091 
— 3.9977 
A^4.6091 
— ©— 5.003 
—0-5.5077 
6.0007 
-♦-6.5009 
—□—7.01 41 
-A- 7.4998 
8.0038 
-5K-8.7928 
-•-9.0056 
9.5132 
10.0031 


x/x t 


Figure 9. Centerline, upper surface, static pressures 
for configuration 1, SPR=0.7. 


The fluidic injection causes an oblique shock just 
upstream of the injection port, as shown in the 
shadowgraph flow visualization in figure 10 and in the 
paint oil-flow visualization in figure 11. Both figures 
are for an NPR of 4.6 with SPR of 0.7. In figure 1 1, the 
location of the shock is identified by the thickened paint 
starting at x/x t = 1.543. In addition, the region of 
separated flow on the flap is evident from the highly 
uneven paint downstream of the injection port, which is 
caused by flow re-circulation. The region of separated 
flow is also indicated in figure 9 by the flat pressure 
profile downstream of the injection port. 



Figure 10. Shadowgraph for configuration 1, 


NPR=4.6, SPR=0.7. 



Figure 11. Paint- oil flow for configuration 1, 
NPR=4.6, SPR=0.7. 

The fluidic injection and resulting oblique shock cause 
an asymmetric pressure distribution as shown in the 
centerline pressures for configuration 1, NPR=4.6, 
SPR=0.7 in figure 12. 



NPR=4.6, SPR=0.7. 



Figure 13. Upper surface, centerline pressures for 
configuration 1 , NPR=4.6. 

Similar trends are witnessed for SPR=0.4 (2% of 
primary mass flow rate) and SPR=1.0 (6% of primary 
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mass flow rate). Increasing SPR increased the strength 
of the oblique shock and moved the shock upstream as 
shown in the upper surface, centerline pressures for 
configuration 1, NPR=4.6 in figure 13. 

Internal flow features are presented for 
configuration 5 (see figure 5) in figures 14 to 17. 
Figure 14 shows the upper surface, centerline, static 
pressures for configuration 5, SPR=0.7 (4% of primary 
mass flow rate). 



♦ 

2.0208 
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2.5072 
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5.005 


5.503 
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6.5098 

-o 

6.9993 

- 

7.5096 


8.0109 


8.7642 


8.9902 


9.4998 

— 

10.0052 


Figure 14. Upper surface centerline pressures for 
configuration 5, SPR=0.7. 

The multiple port fluidic injection causes two oblique 
shocks. The first shock is upstream of the first injection 
port, while the second shock is between the two 
injection ports as shown in the pressure distributions in 
figure 14 and in the shadowgraph for configuration 5, 
NPR=4.6, SPR=0.7 shown in figure 15. 



Figure 15. Shadowgraph for configuration 5, 
NPR=4.6, SPR=0.7. 

It is evident from the flattened pressure profile in figure 
14 that the flow separates at low NPR after the second 
injection port. Flow separation is also evident in figure 
15. The two oblique shocks created by the fluidic 
injection cause an asymmetric pressure distribution as 


shown in the centerline pressures for configuration 5, 
NPR=4.6, SPR=0.7 in figure 16. 



Figure 16. Centerline pressures for configuration 5, 
NPR=4.6, SPR=0.7. 


Similar trends are witnessed for SPR=0.4 (2% of 
primary mass flow rate) and SPR=1.0 (6% of primary 
mass flow rate). Increasing the SPR causes the 
upstream oblique shock to get stronger and move closer 
to the throat and the downstream shock to get stronger 
and move closer to the upstream injection port as 
shown in the upper surface, centerline pressures for 
configuration 5, NPR=4.6 in figure 17. 



Figure 17. Upper surface centerline pressures for 
configurations, NPR=4.6. 

The internal flow features for configurations 2, 3 
and 4 are similar to configuration 5 and were omitted 
for brevity. For further information on the internal flow 
features of these configurations, please refer to 
reference 13. 

Pitch Thrust Vector angle 

Figure 18(a-d) shows a comparison of the pitch 
thrust vector angle for all five configurations at all test 
conditions. 
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8 P (deg) 8 P (deg) 5 P (deg) 5 P (deg) 



NPR 


(a) No Injection (SPR=0.0). 



NPR 


(b) SPR=0.4. 



(c) SPR=0.7. 



NPR 


(d) SPR=1.0. 

Figure 18. Pitch thrust vector angle. 


Translating the upstream injection port closer to the 
throat had varying effects on the pitch thrust vectoring 
capability of the nozzle. In general, increasing the 
number of injection ports from one to two increased the 
pitch thrust vectoring capability at NPR less than 4, 
especially at high SPR. The configuration that 
generated the largest thrust vectoring varied with NPR 
and SPR. In general at SPR =1.0 and NPR<4, 
configuration 5 generates the largest pitch thrust vector 
angle. At SPR=0.7, configuration 4 generally provided 
the largest pitch thrust vector angle for NPR<5. At 
SPR=0.4, all configurations were within 1.5 to 2 
degrees, but configurations 2, 3, and 4 generally had 
higher pitch thrust vector angles than configuration 1 
(single port, baseline injection configuration) at 
NPR<5. At NPR greater than about 4 or 5, there 
appears to be no benefit of the multiple injection port 
concept over the single injection port configuration; for 
some configurations, most notably configuration 5, a 
significant reduction in the pitch thrust vector angle was 
measured. 

The non-linear behavior of the pitch thrust vector 
angle curve is attributed to perturbations caused by the 
injection ports 13 . 

Thrust Performance 

Figure 19(a-d) shows the resultant thrust ratios for 
all five configurations at all test conditions. Translating 
the upstream injection port closer to the throat generally 
has a negligible effect on the resultant thrust ratio 
except at NPR=2 with all SPRs and at NPR<4 with 
SPR=1.0. In the case of the latter exception, 
configurations 4 and 5 increase the resultant thrust 
efficiency. This is an added benefit, since these 
configurations also provided superior thrust vectoring 
over the same NPR range. 

The increased thrust efficiency is caused by a 
shock jump at NPR=4 and is a direct result of the 
separation reduction explained in reference 22. In 
particular, as flow separation is alleviated, the shock 
abruptly moves or "jumps". 

The resultant thrust ratios may be greater than one 
because the ideal thrust is computed based on total 
temperature, total pressure, and weight flow for the 
primary flow only. The measured thrust includes the 
thrust generated from both the primary and the 
secondary flow. The thrust ratio is reported in this 
manner because there was no total temperature 
measurement in the secondary plenum, which is 
essential for computing the ideal thrust of the secondary 
flow. 
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Fr/R Fr/R Fr/R 



NPR 


(a) No Injection (SPR=0.0). 



NPR 


(b) SPR=0.4. 



NPR 


(c) SPR=0.7. 



NPR 

(d) SPR=1.0. 

Figure 19. Resultant thrust ratio. 


Computational Results 

A computational study was conducted on two of 
the experimental configurations: configuration 1 and 
configuration 5 (see figure 5). The configurations were 
run at NPRs of 4.6, 7, 8.78, and 10 all with an SPR of 
0.7. Internal flow features (static pressures and flow 
visualization), pitch thrust vector angles, and thrust 
performance will be presented here in comparison to 
experiment. 

Internal Flow 

Internal flow features are shown for configuration 
1 (single port, baseline injection configuration) in 
figures 20 and 21 and for configuration 5 in figures 22 
and 23. Figure 20(a-d) shows the experimental and 
computational centerline pressures for configuration 1, 
SPR=0.7 (4% of primary mass flow rate). 

The PAB3D-predicted static pressures along the 
upper and lower nozzle surfaces correlated well with 
experimental data at all NPR with a few notable 
exceptions in shock location, strength, and profile. 

First, the shock location on the upper surface at all 
NPRs, and on the lower surface at NPR=4.6, was 
predicted slightly downstream of experimental data. In 
particular, PAB3D predicted the shock at x/xt=1.61 on 
the upper surface compared to a shock location of 
x/xt=1.53 in the experimental data at NPR=4.6 (figure 
20(a)). Also shown in figure 20(a), PAB3D predicted 
the shock at x/xt=1.95 on the lower surface compared 
to a shock location of x/xt=1.89 in the experimental 
data. Second, PAB3D predicted a slightly stronger 
shock on the upper surface compared with experimental 
data for all NPR, as indicated by the larger pressure rise 
in figure 20. On the lower surface, the start of a shock 
is present only at NPR=4.6. Since the shock is not fully 
internal, a comparison of experimental and predicted 
pressure rise is not possible. Lastly, PAB3D predicted 
a steeper pressure profile of the shock on the upper 
surface at all NPR, when compared to the experimental 
data. 

The discrepancies between experimental and 
predicted shock location, strength, and profiles are 
attributed to the thin layer approximation used in the 
computational study. The simplified, uncoupled 
Navier-Stokes approximation would be more 
appropriate to model the complex shock system as 
demonstrated in reference 18. Due to the complexities 
brought about by the secondary stream interactions with 
the supersonic primary stream, a simplified analysis 
utilizing the thin layer approximation was conducted. 
Future studies should investigate the problem utilizing 
the simplified, uncoupled Navier-Stokes approximation. 
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x/x t 


(a) NPR=4.6. 



x/x t 


(b) NPR=7. 



x/x t 


(c) NPR=8.78. 



(d) NPR=10. 

Figure 20. Experimental and computational centerline 
pressures for configuration 1, SPR=0.7. 


Figure 21 shows the predicted flow characteristics 
using CFD Schlieren flow visualization for 
configuration 1, NPR=4.6, SPR=0.7. Qualitatively, 
PAB3D predicted the shock patterns very well when 
compared with the experimental flow patterns shown in 
figure 10. The 6 p =7° thrust vector angle of the primary 
flow is visible in both figures. 



Figure 21. CFD Schileren for configuration 1, 
NPR=4.6, SPR=0.7. 


Figure 22(a-d) shows the experimental and 
computational centerline pressures for configuration 5, 
SPR=0.7 (4% of primary mass flow rate). 

The PAB3D-predicted static pressures along the 
upper and lower nozzle surfaces correlated well with 
experimental data at all NPR with a few notable 
exceptions in shock location, strength, and profile. 

First, the first shock location on the upper surface 
at all NPRs, and the shock location on the lower surface 
at NPR=4.6, was predicted slightly downstream of 
experimental data. In particular, PAB3D predicted the 
first shock at x/x t =1.3 on the upper surface compared to 
a shock location of x/x t =1.24 in the experimental data at 
NPR=4.6 (figure 22(a)). Also shown in figure 22(a), 
PAB3D predicted the shock at x/x t =1.95 on the lower 
surface compared to a shock location of x/x t =1.89 in the 
experimental data. Second, PAB3D predicted a slightly 
stronger first shock on the upper surface compared with 
experimental data for all NPR, as indicated by the 
larger pressure rise in figure 22. On the lower surface, 
the start of the shock is present only at NPR=4.6. Since 
the shock is not full internal, a comparison of 
experimental and predicted pressure rise is not possible. 
Lastly, PAB3D predicted a steeper pressure profile of 
the first shock on the upper surface at all NPR, when 
compared to experimental data. 

The discrepancies between the experimental and 
predicted shock location, strength, and profiles are 
attributed to the thin layer approximation used in the 
computational study as mentioned previously. 
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nd /d nd/d rid/d rid /d 



x/x t 


(a) NPR=4.6. 



x/x t 


(b) NPR=7. 



Figure 23. CFD Schileren for configuration 5, 
NPR=4.6, SPR=0.7. 


Figure 23 shows the predicted flow characteristics 
using CFD Schlieren flow visualization for 
configuration 5, NPR=4.6, SPR=0.7. Qualitatively, 
PAB3D predicted the shock patterns very well when 
compared with the experimental flow patterns shown in 
figure 15. The 6 p =6.5° thrust vector angle of the 
primary flow is visible in both figures. 



x/x t 


(c) NPR=8.78. 



(d) NPR=10. 

Figure 22. Experimental and computational centerline 
pressures for configuration 5, SPR=0.7. 


Pitch Thrust Vector angle 


Figures 24 and 25 show the experimental and 
computational pitch thrust vector angle for 
configurations 1 and 5, respectively. 



Figure 24. Experimental and computational pitch thrust 
vector angle for configuration 1, SPR=0.7. 


The PAB3D-predicted pitch thrust vector angles 
were within the error band of the experimental data as 
shown in figures 24 and 25. 
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Figure 25. Pitch thrust vector angle for configuration 
5, SPR=0.7. 


Thrust Performance 

Figures 26 and 27 show the experimental and 
computational resultant thrust ratios for configurations 
1 and 5, respectively. PAB3D predicted the thrust ratio 
within 0.57 and 0.37 percent of the experimental thrust 
ratio for configurations 1 and 5, respectively. The 
predicted thrust ratios are within the experimental error 
band as shown in figures 26 and 27. 



Figure 26. Resultant thrust ratio for configuration 1, 
SPR=0.7. 



Figure 27. Resultant thrust ratio for configuration 5, 
SPR=0.7. 


Conclusion 

The objective of the study was to investigate 
multiple injection ports in a two-dimensional, 
convergent-divergent nozzle for fluidic thrust vectoring 
through experimental and computational means. In 
particular, a two-dimensional, convergent-divergent 
nozzle with two injection ports was hypothesized to 
enhance the thrust vectoring ability of a two- 
dimensional, convergent-divergent nozzle with a single 
injection port without increasing the secondary flow 
requirements and incurring significant performance 
penalties. In addition to verifying this hypothesis, it 
was the goal of the study to validate the computational 
fluid dynamics code, PAB3D, for multiple injection 
ports in a two-dimensional, convergent-divergent 
nozzle. The validation of the computational fluid 
dynamics code would allow future investigations to be 
performed solely by computational means. 

The study included testing five experimental 
configurations and running computational simulations 
on two of the experimental configurations. The 
experimental configurations were run at nozzle pressure 
ratios up to 10 with secondary pressure ratios up to 1. 
The computational configurations were run at nozzle 
pressure ratios up to 10 with a secondary pressure ratio 
of 0.7. 

Results of the experimental study indicate that 
multiple injection ports increase thrust vectoring with 
increased performance benefits at nozzle pressure ratios 
less than 4 with high secondary pressure ratios, 
compared to one injection port. At nozzle pressure 
ratios greater than 4, there appears to be no benefit of 
the multiple injection port over one injection port. To 
take advantage of improved thrust vectoring without 
increasing flow requirements, a nozzle could be 
designed with location-controllable injection ports that 
would optimize vectoring over a range of nozzle 
pressure ratios. 

Results of the computational study indicate 
PAB3D is valid for predicting the internal performance 
quantities of a two-dimensional, convergent-divergent 
nozzle with multiple injection ports. For future studies, 
PAB3D could be used with confidence in exploring 
different configurations for multiple injection ports. 

This study investigated a very limited number of 
configurations to verify the hypothesis and no attempt 
was made to optimize the configuration. Future studies 
are needed to obtain an optimal configuration for 
various NPR conditions. PAB3D could be used in 
conjunction with a computational optimizer to find an 
optimal configuration. 
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